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Using the known connection between Schroedinger-like equations and Dirac- 
hke equations in the supersymmetric context, we discuss an extension of FRW 
barotropic cosmologies in which a Dirac mass-like parameter is introduced. New 
Hubble cosmological parameters Hk(t?) depending on the Dirac-like parameter are 
plotted and compared with the standard Hubble case Ho (77). The new Hk(77) are 
complex quantities. The imaginary part is a supersymmetric way of introducing 
dissipation and instabilities in the barotropic FRW hydrodynamics. 



1 - Introduction 

The scale factor a(t) of a FRW metric is a function of the comoving time t 
obeying the Einstein-Friedmann dynamical equations supplemented by the 
(barotropic) equation of state of the cosmological fluid 



(P + 3p) , (1) 



a 3 

_ SttGp _ K 
"3 a2 ' 

P = (7 - (3) 

where p and p are the energy density and the pressure, respectively, of the 
perfect fluid of which a classical universe is usually assumed to be made 
of, K = 0,±1 is the curvature index of the flat, closed, open universe, re- 
spectively, and 7 is the constant adiabatic index of the cosmological fluid. 
A simple Riccati route of solving the system of eqs. (l)-(3) proposed by 
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FaraoniP has been used by Rosu to develop a factorization (nonrelativistic 
super symmetry) approach of the barotropic FRW cosmologies)^ in which a 
simple explanation for a currently accelerating universep^ in terms of 'Dar- 
boux fluids' with effective adiabatic indices is possible. 

In this work, we first briefly review the supersymmetric factorization 
methods for barotropic FRW cosmologies in section 2. Next, in section 3, 
we present Dirac-like (first-order) coupled differential equations and their 
associated second-order differential equations. This allows a simple exten- 
sion of barotropic FRW cosmologies that include a Dirac mass-like constant 
parameter K. 

2 - Supersymmetric methods 

Combining the equations (l)-(3) and using the conformal time variable r] 
defined by dt = a(r])dr] one gets the equation 

= . (4) 

where c = |7 — 1. The case k = is directly integrable,^ and will be skipped 

henceforth. The logarithmic derivative Hq = ^ is Hubble parameter of 
the universe which is a fundamental cosmological quantity ever since Bub- 
ble's discovery of the expansion of the universe inl929. Interestingly, passing 
to the function Hq in Eq. the following Riccati equation is obtained 

Ho cHg Kc = . (5) 
Employing now Hq = one gets the very simple second order differential 
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W^ -1- KC Wk = . (6) 

For K = 1 the solution of the latter is wi = Wi cos(cr/ -|- d) , where d is 
an arbitrary phase, implying ai{ri) = Ai[cos(c77 + d)]^/'^ , whereas for k = 
— 1 one gets w_i = W_isinh(c7/) and therefore a_i(r/) = A_i[sinh(cr/)]^/'^ , 
where W±i and A±i are amplitude parameters. 

The point now is that the particular Riccati solutions H(j" = — tan crj and 
Hq = cothcr/ for k = ±1, respectively, are closely related to the common 
factorizations of the Schrodinger equation. Indeed, Eq. (jH)) can be written 

w — c(— Kc)w = (7) 

and also in factorized form using Eq. (6) one gets 

+ cHo) (-^ - cRo) w = w" - c(Ho + cHg)w = . (8) 
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To fix the ideas, we shall call Eq. (jH)) the bosonic equation. On the other 
hand, the supersymmetric partner (or fermionic) equation of Eq. ((Hj) will be 

- cHo^ + cHo^ Wf = Wf - c(-Ho + cH§)Wf = wj-' - C • CK,susyWf = 

(9) 

Thus, one can write 

2 _ J c(l + 2tan^cr/) if k = 1 



C«;,susy(??) = -Hq + cHf, 



c(- 1 + 2coth^ cr?) iffi: = -l 



for the supersymmetric partner adiabatic index. The solutions are Wf = 
— , , ,^ and Wf = ^-r^^ — r for K = 1 and k = — 1, respectively. 

3 - Dirac-like formalism 

The Dirac equation in the susy nonrelativistic formalism has been discussed 
by Cooper et al^ already in 1988. They showed that the Dirac equation 
with a Lorentz scalar potential is associated with a susy pair of Schroedinger 
Hamiltonians. This result has been used later by many other authors. Here 
we make an application to barotropic FRW cosmologies that we find not to 
be a trivial exercise except for the uncoupled 'zero-mass' case (subsection 
3.1). 

3.1- Let's introduce now the following two Pauli matrices 

a = -iay = -i(^^. and /3 = = 1^ ° J ^ 

and write a cosmological Dirac equation 

T^FRW^ ^ ^ ax(icHo)]W = , (10) 

where = ( | is a two component 'zero-mass' spinor. This is equiva- 

V ^2 J 

lent to the following decoupled equations 

Dr^Wi + cRqWi = (11) 
-D^W2 + cHot«2 = . (12) 

Solving these equations one gets wi oc l/cos(cr/) and ■W2 oc cos(c?7) for k, = 1 
cosmologies and wi oc l/sinh(cT/) and W2 oc sinh(c77) for k = —1 cosmologies. 
Thus, we obtain 

\ W2 \ Wb / 
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This shows that the matrix 'zero-mass' Dirac equation is equivalent to the 
two Schroedinger equations for the bosonic and fermionic components. 



3.2- Consider now a "massive" Dirac equation 

H^RW^ ^ [^^D^ ^ cJx(icHp + K)]W = KW 



(13) 



where K may be considered the mass parameter of the Dirac spinor. Eq. 
is equivalent to the following system of coupled equations 



D^wi + (icHo + K)wi = Kw2 
-D^W2 + (icHo + K)w2 = Kwi . 



(14) 
(15) 



These two coupled first-order equations are equivalent to second order dif- 
ferential equations for the two spinor components: 
The fermionic spinor component can be found directly as solutions of 



and 



D^wr 



c^ ( 1 + 2 tan^ cij) + 2icK tan ct] 



c^(— 1 + 2coth^c?7) — 2icKcothc?7 
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for K = 1 



for K 



whereas the bosonic components are solutions of 



D^W2+ + 



and 



D^w^ + 



2icK tan cr/ 



c + 2icKcoth CT] 
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for K = 1 







for K. 
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The solutions of the bosonic equations are expressed in terms of the Gauss 
hypergeometric functions 2F1 in the variables y = e''^^ and y = e*^^, respec- 
tively 



w+ = Ay p 2F1 [ - ^(p + iq); -^(p - iq), 1 - p; -y^ 

ByP 2F1 [^(p - iq), ^(p + iq), 1 + p; -y^ 



+ 



(20) 



and 



W2- = C(-l)-2'^y-i'-2Fi 



^(r + s),-^(r-s),l 



r;y 



+ 



D(-l)i'-ySFi[^(r-s),^(r + s),l + r;y2' 
respectively. The parameters are the following: 



(21) 
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whereas A, B, C, D are superposition constants. 

Based on these K- modes, we can introduce bosonic Hubble parameters 
depending on the parameter K 

and 

and similarly for the fermionic components by changing w^ to w^ in eqs. (|22|1 
and ((221), respectively. 

In the small K limit, K — > 0, the ordinary FRW barotropic cosmologies 
are obtained. The Hubble parameters corresponding to the K- dependent 
bosonic modes are plotted in the Figures (1) -(4). 

4 - Conclusions 

We come now to the interpretation of the mathematical results that we 
displayed in the previous sections. The parameter K introduces an imagi- 
nary part in the cosmological Hubble parameter H. Since the latter is the 
logarithmic derivative of the scale factor of the universe one comes to the 
conclusion that the supersymmetric techniques presented here are a way to 
consider dissipation and instabilities of barotropic FRW cosmologies. 
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Figure 1: The real part of HJ^ obtained from the logderivative of the bosonic 
mode w^(y; i, i) for rj £ [0, 10] and K G [0, 10] in the case of a radiation 
dominated universe (c = 1) . 




Figure 2: The imaginary part of Hj^ for the same case. 
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Figure 4: The imaginary part of Hj^ for the same case. 



